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We derive a criterion to determine when a translationally invariant matrix product state (MPS) 
has long-range localizable entanglement, where that quantity remains finite in the thermodynamic 
limit. We give examples fulfilling this criterion and eventually use it to obtain all such MPS with 
bond dimension 2 and 3. 
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Localizable entanglement (LE) PQ is a multipartite 
measure characterizing the (maximal averaged) bipartite 
entanglement present in a system. Motivated by quan- 
tum repeaters [5] , where it provides the natural figure of 
merit, the LE has also been applied to many-body physics 
problems [3HZ], whereby it can reveal hidden correlations, 
not detectable by standard observables [3J. 

In both scenarios, there is an underlying spatial struc- 
ture. Hence, peculiarities are expected whenever a finite 
amount of entanglement between two particles with arbi- 
trary distance can be created by properly measuring the 
rest of the particles. Such states are said to have long- 
range localizable entanglement (LRLE) [8] and play an 
important role both in the context of quantum repeaters 
and spin chains. In the first one, they are those for which 
entanglement can be established at arbitrary distances. 
In the latter, phase transitions are signaled by a finite 
value of the LE between any two particles irrespective of 
their distance. States with LRLE thus play a crucial role 
in those two contexts. But, which are those states? 

As most multipartite entanglement measures, the LE 
is very hard to determine for general states. A notable 
exception is the set of matrix product states (MPS) [S]. 
This family of states describes the ground state of gapped 
ID spin chains as well as those states created by sequen- 
tial generation, as it is the case of atoms in a cavity [TU] , 
They are characterized by a set of rank-three tensors, 
each one associated to a spin. One of the indices cor- 
responds to the spin in the z-basis, and the other two 
run from 1 to D, where D is called the bond dimension. 
MPS of arbitrary bond dimension arc dense in the set of 
all multipartite states |11) . and thus, they are very rel- 
evant to describe many-body systems. For translation- 
ally invariant systems and low bond dimensions (D = 2) 
the LE can be determined exactly (3J. For larger bond 
dimensions, one may find relatively tight lower bounds 
using Monte-Carlo Methods [T21H5] . 

In this paper we fully characterize translationally in- 
variant MPS with LRLE for arbitrary bond dimension. 
In particular, we give a set of necessary and sufficient 
conditions for such states. As we show, those conditions 
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can be turned into a set of polynomial equations, and 
thus provide us with a precise criterion to determine if a 
state has LRLE or not. Furthermore, we give examples 
of non-trivial states with that property and provide the 
full sets of MPS with LRLE for bond dimension D = 2 
and D = 3. 

The LE is defined as the maximum average entangle- 
ment that can be generated between two spins of a spin 
chain by measuring the remaining ones [1] . Let p denote 
the density matrix of the original state. With probability 
Pi the outcome of a measurement M. will be i and the 
system will be in the corresponding two-particle state p\. 
Hence, the LE is given by 

L c > E (p) = sup ^p^(^), (1) 
Mec j 

where C is the class of allowed measurements and E(-) 
an entanglement measure. 

Our system of consideration is an open chain of N spin- 
s' particles along with two auxiliary particles of spin S' 
at each of the boarders. The N particles of the actual 
chain are the ones to be measured (measurement out- 
comes i = (ii, ijv))j an d the class of allowed mea- 
surements C is the set of local projective von Neumann 
measurements, where the same measurement is carried 
out on each party (in particular, we exclude adaptive 
strategies). Therefore, the maximization of the average 
entanglement is performed by choosing the optimal phys- 
ical basis Ui)} 2 ^ 1 - 

The question to be answered in this letter is for which 
translationally invariant MPS a finite amount of entan- 
glement can be localized between the two ancillas in the 
limit N — > oo. We assume the state of the system to be 
translationally invariant apart from boundary effects; for 
this reason the rank-three tensors corresponding to the 
spin-S particles are taken equal. Those consist of com- 
plex D x D matrices Aj (i = 1, d = 25+1) and 
can be taken to be in canonical form, in which the maps 

£P0 = Eti A i XA l and 8(X) = ELi4* A i satisfy 
(cf. i) 

£(I) = I, £(A) = A (2) 

for some diagonal positive definite matrix A. Our goal is 
to find necessary and sufficient conditions on the matrices 
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FIG. 1. Spin chain with N spin-S 1 particles and one auxiliary 
particle of spin 5' at each of the borders. The matrices of the 
real particles are A, and those of the auxiliary particles 
and Q, respectively. 



{Ai}f =1 to give rise to LRLE for some matrices of the 
auxiliary particles. Those can be chosen at will and are 
denoted by P,Q : C D ' -> C D , where D' = 25' + 1 < D 
is the Hilbert space dimension of the individual auxiliary 
spins. The initial MPS is therefore 

D' d 

|V>)=E E (k\^A il ...A iN Q\l)\i 1 ...i N )®\k,l), 



k,l=l ii ,. ■ -,iN — 1 



(3) 



where the Hilbert space vectors of the auxiliary par- 
ticles are denoted by round brackets, c.f. Fig. [T] 
Subsequent to a measurement Ai the initial state of 
the system p — reduces with probability pi — 

tr(^A ll ...A lN QQ^Al N ...Al i V) to pf 4 = h^X^I, 
where (excluding cases with pi = 0) 
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(4) 



k,l=l 



is the normalized two-particle state after the measure- 
ment. In Appendix [A] it is shown that any MPS with 
LRLE for D' > 2 also has LRLE for D' = 2 (the con- 
verse is obvious), and also that w.l.o.g. we can take P 
and Q as isometries. We can thus choose D' = 2 and 
the concurrence [14] as the measure of entanglement, 
E(pf A ) = 2|det(*i)|, where <Fi is the matrix with co- 
efficients = ^(fclPt^.-.A^QIO- Dropping the 
superscripts C and E, the LE reads 



L(P) 



2 sup 

{!*>} 



E 



detfPUv, ...Ai 



(5) 



This optimization problem is in general hard, since the 
sum needs to be evaluated for large N before being able 
to optimize over the physical basis {|«)}f =1 . Its maxi- 
mization succeeded only in special cases, like for a chain 
of N spins without auxiliary ones and measurement on 
all spins but those at the borders [H]. In this case 
Tl^ 1 Ei-=i I det(A ij .)| 2 / D = 1, which implies that there 
is LRLE for Ai = onUi. In the following we slightly gen- 
eralize this family of states and adapt them to our system 
that includes the auxiliary spins at the boundary: 
Example 1. — Block structure of unitaries. 



(6) 



where Pi is a q x q permutation matrix. The Uf are 
n x n unitaries (D — nq) and <x? > V k = 1, q, i = 
1, d. We can realize that this MPS has LRLE for 
n = 2 by choosing pt = | f)(l| + | \){2\ (where | f), I i) 
are basis vectors of the auxiliary qubits and |1), |2), ... 
the basis vectors of the basis in which the matrices 
are given) and Q = D~W ^f=i (U)(t I + (-l) J '|j)(4. I)- 

Then, due to 8(1) — I and therefore £? =1 \ a i\ 2 = 1 V k = 
1, D/2 ([5| is finite in the thermodynamic limit. The 
fact that there is also LRLE for n > 2 will become clear 
below. 

In the above examples the matrices Ai exhibit a block 
structure of unitaries. One may think that all states with 
LRLE need to have this property in some basis {|i)}f =1 - 
Interestingly, this is not the case as shown by the follow- 
ing counterexample for D = 3: 

Example 2. — Non-unitary matrices. 



A x = 



1 

1 

1 



A, = 




A,=-. 



1 

1 
1 



(7) 



with Pt = 



t)(l| + U)(2|andQ = i[|l) + |2) + |3)] (t | 
+ [|2) — |3)] (I |. It is a simple exercise to show that 
(|5) remains finite in the limit N — > oo. The matrices Ai 
fulfill £ (I) = I. However, £ (I) ^ I, and thus they are not 
of the form Ai — ctiUi in any basis {|«)}f =1 . Note that the 
MPS can numerically be verified to be injective [9], and 
thus it is the unique ground state of a local translationally 
invariant gapped Hamiltonian. Moreover, it can be easily 
shown (analytically) that it is invariant under the Z<i- 
symmetry generated by the transformation |1) — » |1), 

|2>->|2),|3)->-|3>. 

The question arises of how one can check whether a 

given MPS has LRLE without having to resort to eval- 
uating ([5| for large N numerically. In the following a 
necessary and sufficient criterion will be derived, which 
allows to decide this based on the matrices {Ai}f =1 di- 
rectly. We first rewrite ^ by inserting to the right of 
Ai . a projector on the subspace spanned by the row vec- 

. P jt ■ , 



tors of F' A il ...A ij , which can be written as 

where , : C 2 -> 
now 



is an isometry. Thus, ^ reads 



L(p) = 2 sup £ |det(PtA ll P J 1 i )| ]T |det(P l 1 U 42 P 2 i!l2 



my 

E 



det(P£ 



1 lN-l Ai » <! 



k=l 



) 

(8) 

Using the SVD of P^ 11 ,- Aj.Pf - , which is 

diag(cr^ 1) , cr^ 2) ) Vl (neglecting the indices it, . . . , ij-x), 
along with the inequality of arithmetic and geometric 



means, one obtains for the factors of (|8 

y|det(pr lt , a-.p? ,) 

- 2 X! tr i^il^li j -l A 'h V li,...,i^ 3 ii,...,i j A \ ] V li^..,i ] -i 



(9) 



since J^i A ij A l = The inequality becomes an equal- 

ity if and only if = a ^ , i.e., H^-i^K,.,^ is 
proportional to a unitary. The intuitive fact that all fac- 
tors but the last of ^ have to be exactly 1 to get LRLE 
is proven in Appendix [B] Therefore, a necessary condi- 
tion for any MPS to give rise to LRLE is that there exists 
an isometry Y\ such that for a certain basis 
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A: 



oc U ll ... i] V ix, 



., d, 



(10) 



denoting some 2x2 unitary. Redefining 
Pi .-. — > Pi i Vi 1 ...i i shows that one can require 



the RHS of (10 1 to be the identity. After multiplying 
this from the right by P^ £ . , one obtains 



j-if 



A 



<'t 



which can also be written as 



(11) 



P f A,...A- OcPf f 

t l , 



(12) 



As illustrated by Example 1 and 2, Q can always been 
chosen such that ^ is finite, which is shown rigorously 
in Appendix [Cj Criterion (12) is thus necessary and suf- 
ficient for the emergence of LRLE [TB] . 

If we define Pt := | \){x\ + \ l)(y\,i.e., (x\x) = (y\y) = 1 



and {y\x) = 0, (12) is equivalent to 



{x\A ll ...A l] A\ j ...A\ i \x) = (y\A ix ...A ij A\....A\\y), (13) 
(y\A il ...A ii 4....Al\x)=0. (14) 

After further defining V := \x)(x\ — W ■= \x)(y\, 

the map Si(X) := A\XAi and 



V,, 



:= A\....A\\x){x\A il ...A ij 

= £ lj o...o£ H (V), 

:= A\ j ...A\\x){y\A il ...A i . 

= I i} o...o£ n (W), 



we see that ( 13 ) and (H 4| are equivalent to 



tv(V iu ... <ij )=tv(£ ij o...oS il (V))=0, 
ti{W il ,... tij )=ti(£ i .o...o£ il (W)) = 0, 



(15) 
(16) 
(17) 
(18) 

(19) 
(20) 



respectively. Thus, if we define the subspace S 
span{V, W, Vi u ., 



-\i :• @ and 

( |20| indicate that the occurrence of LRLE is equivalent 
to tr(S*) = V £ G S. By definition S is closed under the 
application of any which leads to the following cri- 
terion characterizing translationally invariant MPS with 
LRLE: 

Theorem 1. — For the MPS considered here, there is 
LRLE if and only if for a certain basis {|«)}f = i there 
exists a subspace S of the vector space of D x D matrices 
satisfying the following conditions 

1. S is closed under all linear maps £\ defined as 



£i(X) = A}XAi, i.e., £ t {S) C S V * = 1, 



d. 



\x){x\ - \y){y\, W 



2. tr(S) = V S e S, 

3. 3 |a;), \y) £ C D s.t. V 
\x)(y\eS. 

Note that if Theorem 1 is fulfilled, £ h o...o£ h (W^) = 

W\ t j. is also traceless, i.e., S could additionally be 

required to contain W* = \y)(x\ and correspondingly to 
be equal to its adjoint, V S e S => S 1 ^ e S. 

Theorem 1 can be used numerically to determine, 
whether a given MPS has LRLE, since it imposes con- 
ditions on the matrices {A}f=iJ which can be repre- 
sented by a set of polynomial equations: The entries of 
{A}f = i will give rise to the coefficients of those equa- 
tions, whereas all other quantities introduced below will 
constitute their variables to be determined numerically. 
The first variables to be introduced are orthonormal basis 
vectors {S k } 7 ^ =1 of S , where n < D 2 — 1 is the dimension 
of S. Condition i is equivalent to requiring that any basis 
vector S of S is mapped by any £ i into S, i.e., for all 
k = 1, n 



I. A\ S k Ai = a>i' l S l V i = 1, d, 



i=i 



tv(S k ^S l ) = 4.z V I = 1, 



(21) 
(22) 



where {a 4 ' } i= L 1 l '"^ n are scalar complex variables of the 

set of equations to be solved numerically. The {A}f =1 
are the matrices in a possibly different physical basis 
{|*)}i=i- h can t> e stated as 



II. tr{S k ) = V k = 1, 



(23) 



Furthermore, since according to condition iii V = \x)(x\ 
\y)(y\ and W = \x)(y\ have to be also contained in S, 



iii. J> fc s fc H*)(*|- \v)(y\, 

k=l 

n 

"£w k S k = \x)(y\ 



(24) 
(25) 



k=l 



with {v k }^ =1 , {w k }^ =1 , and the coefficients of \x) and 
\y) as other scalar complex variables. Last, a rotation in 
the basis of measurement is implemented by 



IV. 



A i = ^U ij A j , UU^ = 



-dxd> 



(26) 



constituting the last of the set of equations (21 1 to (26) 
to be solved. Generally, a set Pi{x\, x m ) = of s 
polynomial equations (i = 1, s) with m variables 
{xj}JL 1 can be solved by means of a Grobner basis [T7] , 
A Grobner basis {<7t(xi, x m )}f =1 is a special basis 
in the vector space of functions involving the variables 
{xj}'fL 1 : It has the property that the set of equations 
to be solved is equivalent to the set gi{x\, x rn ) = 
(i = 1, s), which, in contrast, can be solved by back- 
substitution while having to deal with the solution of 
polynomials involving only one variable at a time. For 
instance, one of the new equations might involve only 
X\. After solving it numerically the result can be in- 
serted into another equation involving, e.g., only x\ and 
Xi etc. (cf. Gaussian elimination). A Grobner basis can 
be found systematically by use of Buchberger's algorithm 
[T7] , which is doubly exponential in the complexity of the 
set of equations to be solved. In our case, this implies 
a computational cost that is doubly exponential in the 
square of the bond dimension. However, D is a constant 
and in particular independent of the length of the spin 
chain. If a simultaneous solution to ( 21 ) to ( 26 1 is found, 
the MPS has LRLE. If even for n = D 2 — 1 no solution 
is found, it does not. 

We now employ Theorem 1 analytically to determine 
the complete sets of MPS with LRLE for D = 2 and 
D = 3. In the former case, we will reproduce the finding 
that the matrices need to be proportional to unitaries. In 
the latter, we obtain the result that matrices of the type 
of Example 2 are the only non-trivial matrices which give 
rise to LRLE. 

In both cases we take S to be equal to its adjoint, 
i.e., it has to contain the matrices V = \x)(x\ — \y)(y\, 
W = \x)(y\ and WL For D — 2 this implies that S is 
the full subspace of traceless 2x2 matrices. If we take 
any element S e S, ti(AtSAi) = shows that A { A\ 
is orthogonal to S with respect to the Hilbert-Schmidt 
scalar product. However, the orthogonal subspace of S 
is spanned by the identity, i.e., AiA\ oc I for any MPS 
with LRLE. 

For D = 3 we consider first the case of S = 
span{V, W, W^}. We define r = span{|a;), \y)}, whereas 
tracelessness of W and V implies (x\y) — and (x\x) — 
(y\y), respectively (the latter will be taken equal to 1 in 
the following). It follows that S C r x r. If w.l.o.g. we 
set |x) = |1) and \y) — |2), we observe that we retrieve 
the case of D — 2, as 



.4 



15=3 



A°= 2 
B, c 4 , 



where Af =2 is proportional to a unitary, and Bi £ C lx2 
and Cj G C are arbitrary. If S is of larger dimen- 
sion, there are other matrices V — \x')(x'\ — \y')(y'\, 
W = \x')(y'\ and (e.g., V = £ t (V) with \x') = 
A\\x), \y') = A\\y) etc.), which do not span the same 
space as V, W and . We define r' = span{|cc'), \y')} 
and take \x') and \y') as orthonormal (which is possi- 
ble due to tr(V') = tr(W') = 0). Since t and r' inter- 
sect, we can assume the intersection to be w.l.o.g. along 
\y). As y', W' and W 1 * span the full space of trace- 
less matrices contained in r' x r', we are allowed to ro- 
tate the basis {\x'), \y')} for r' such that \y) — \y'). 
Then spanjF, W, W^} and span{V, W, W'^} might 
only intersect in |y)(y|, which is, however, not trace- 
less. Consequently, S is at least of dimension 6. In the 
case of S = span{F, W, W\ V, W, W'*} it is simple 
to construct the orthogonal complement of S, which is 
S = span{I, |n)(n'|, |n')(n|}, where \ri) and |n') are the 
normal vectors of r and r', respectively. This structure 
of S uniquely defines r and r' and vice versa. There- 
fore, another t" ^ r, r' would not be consistent with 
this S, and the only possibility would be S — span{I} 
corresponding to AiA\ oc I (i.e., dim(<S) = 8). Thus, 
dim(5) = 6 is the only remaining case which might lead 
to a non-trivial MPS. In this case there exist the two 
subspaces t and r' of the afore-mentioned type. Based 
on their properties an ansatz for Ai is made in Appendix 
iDl which after an elaborate case differentiation leads to 



Ai oc [e*|l) + e^|3)](Z|+e^'|2)(m|, (28) 



with I, m = 1, 2, 3, I ^ m, as in Example 2. 

Within the framework of MPS we have specified the 
localizable entanglement of a spin chain with one auxil- 
iary spin at each of the borders. We have shown that 
LRLE can be detected by placing qubits at the ends of 
the chain. Based on that we were able to derive a theo- 
rem according to which it can be checked directly from 
the matrices of the MPS, whether it possesses LRLE. 
How this can been done in practice has been indicated 



by eqs. ( 21 ) - ( 26 ) , which is a polynomial set of equations 



(27) 



that can be solved numerically in a systematic manner. 
Furthermore, we provided non-trivial examples of MPS, 
for which those equations have a simultaneous solution, 
determining the full sets of MPS with LRLE for D = 2 
and D = 3. 
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Appendix A: Proof that LRLE for D' > 2 implies 
LRLE for D' = 2 

We would like to show that if the LE acquires a finite 
value for some D' > 2, isometries P, Q G P := {P : 
C 2 C 13 s.t. Ftp = I 2x2 } can be chosen such that ^ 
is non-vanishing in the thermodynamic limit. 

We define the coefficient matrix of IV'i ) via 
(*l)w = ^(fclP^n-^QIO (here P, Q : C D ' -» C c 

Let Pi and 



can be arbitrarily chosen) 



denote the 



isometries contained in V projecting A il ...A iN into the 
subspaces of dimension 2 corresponding to the maximum 
Schmidt coefficients. For the resulting state we obtain 
the coefficient matrix $j = -^P|'v4j I ...Aj Jv .Qi, where pi 
is the probability which arises from keeping track only of 
the projection into the space of the two highest Schmidt 



coefficients, i.e., pi = tr(P 
Now, we want to show that 



Ai. ...Ai. 



..4 



A lower bound for the RHS of ( |A3[ ) can be found by us- 
ing the entropy of entanglement of the qubit system into 
which has been projected, E(^i) = J"(C(#i)) < C(#i) 
[H], where P(C) = h((l + s[\ - C 2 )/2) and h(x) = 
— xlog 2 x — (1 — x) log 2 (l — x), and thus P(C) < C for 
< C < 1 . In our case by definition of Pi and Qi the 
entropy of entanglement is 



> 



Ef^A = ^ — fog, — - — 

V ; Pi + P 2 &2 Pi+P 2 

Pi(D'-l) 
P X {D> -l) + l-P 1 

I -Pi 
Pi(D> - 1) + I 
r(PuD') 



Pi . Pi 
P1+P2 ° g2 Pi + 



+ P-2 

Pi( D'-l) 

1) + 1 - P 
I -Pi 



Fi l0g2 Pi(Z?'-f) + I-P 1 



(A5) 



since for fixed P\ the uppermost sum is minimized for 



P 



_ 1-Pi 



Therefore, we have 



(Al) 



where E(-) is the entropy of entanglement [18] . C(-) the 
concurrence [14] and f(D') some finite- valued function to 
be calculated below. E(-) is given by the von Neumann 
entropy of the normalized reduced density operator of 
the bipartite state pab, E(pab) = — ^(pa log^P/O) with 
Pa — t^sipAB)- The concurrence, C(-), is defined for 
a pure state \ip) of two qubits as C(\ip)) — \(ifj*\a y ® 
cr y \i/})\ = 2\ det(^)| < I, where * is the 2 x 2 coefficient 
matrix of the state. 

We drop the index i for the variables introduced in 
the following and denote by Qi > Q 2 > ... > Qn' the 
eigenvalues of the unnormalized 1-particle reduced den- 
sity operator, p' A , of p' AB = Pi\i^){^\ = Pi Pab after 
the measurement. Thus, we get 



£7(*i) < w(P u D') 



C(*i) " r(P u D<) 



(A6) 



which is finite for < P% < 1- In the limit P x — > 1 
the ratio takes the value D' — 1, thus verifying that 
there exists a finite-valued function g(D') such that 
£(*i)/C(*i) < g (m. 

Now we apply ((AlJ) to bound the LE for arbitrary D' 

L(p) = sup VpiP^O < f(D') sup VpiC(*i) (A7) 

with ^ ; = ^JPt A^ ...A iN Qi. In analogy to ^ we obtain 

L(p) < 2f(D') sup |det(P|A il ...A ijv Q 1 )|. (A8) 

{!''>} i, t„ 



D' D , 

k=l 



and it is sufficient to show that 



D' 



Em < g(D')C(h). 



(A2) 



(A3) 



In order to do so we first derive an upper bound for 
the LHS of ( |A3| . The eigenvalues of the normal- 
ized 1-particle reduced density operator, pa, of pab — 
l^i M )(^i M | are P k = Q fe /£ n Q„. Hence, we obtain 



D' 



E(Vi) = E(pab) = -J2 Pk lo & Pk 



fe=i 



< -Fx log 2 P x - (1 - P x ) log 2 
:=w(Pi,D>), 



(A4) 



since the sum from the second term on is maximized for 
P 2 = - = Pd- = (1 - Pi)/(D' - 1). 



We introduce an e-net (e > 0) A/" e of isometries contained 
in V 1 that is, for any such isometry P there exists a V G 
M e with operator norm ||P — V|| < e. By choosing e 
sufficiently small we get up to orders in e 

L(p) < 2/(P') E E E |det(VU n ...A 1JV W)| 

x ^V.Pi^W.Qi (A9) 
:=2/(Z>') £ ]T w(V,W). (A10) 



If the LE is non-zero in the thermodynamic limit, 
(A10) is lower bounded by L(p) > 0. Consequently, 
for all r ^ e there must exist r-regions IZl, 7£ 2 C J\f e 
containing isometries V, W such that V V £ IZl and 
VWeR^we have ||V - V'|| < r and ||W - W'|| < r, 
respectively, with 



V{D') E E ^V,W)>L(p)^^, (All) 



V6RJ WgRj 



where Vr. denotes the volume, i.e., the number of points 
in the e-discretization of region 1Z. By choosing r small 
enough, one obtains w(V ,W) = lj(Y , W')(l+C(r)), such 
that sum in the LHS of (All I contains Vri V^2 /V?f con- 
stant terms, and therefore 



2/(D>(¥',W')(l + 0(r))>L(/)). 



(A12) 



Hence, for sufficiently small r one gets w(V,W) > 0, 
i.e., ([5| is non- vanishing for some isometries V and W, 
if there is LRLE for some D' > 2 □ 



Appendix C: Proof that (121 is also a sufficient 



condition 

Here we show that Q can always be chosen such that 
^ is finite whenever ( fl2| is fulfilled. We denote the 
proportionality factor in ( |12[ ) by 7^ for which £(T) = 

Eti A i A \ = 1 implies £\ \liu-,ij\ 
thus reads 



L(p) = 2 J2 l^K 



JVt 



2 = h u ...,^A 2 - (§ 
lii,..^! 2 - (ci) 



Appendix B: Proof that all factors but the last one 
of |8| have to be 1 

We want to demonstrate that in the limit N — > 00 the 
sum ([8| can be non-zero only if one can choose such 

that for all well-defined P; - j (i.e. those for which 

i \ ? ■ ■ ■ 1 l j — 1 x 

F J< A ll ...A i ._ 1 ^ 0) the inequality 



< 1 



(Bl) 



is an equality. 

Formally this claim reads 

LRLE =^ 3P f s.t. V i lt = 1, d 

E|^(<:.^-i^X,..,0| = 1 . ( B2 ) 

or equivalently, 

-1 LRLE 4= VP 1 3 ij_x) s.t. 

E| det ( p ^- 1 ^- 1 ^ p -.-.0k 1 ' (B3) 

which is what we want to show in the following. Thus, 
we assume that for the sum 

E |det(P t A il ...A i .Pf 1 := 1- A^(P f ) (B4) 

i\ ,...,z s 

there is a minimum integer such that 

minA^Vpt) := A, > 0. (B5) 

pt 

After building blocks of s terms in ^ , each of them can 
be upper bounded by 



l-A(P{t )<1-A, 



which shows that the LE must fulfill 

L(p)<(l-A^- 1 ^, 0, 



i.e., there is no LRLE. 



(B6) 

(B7) 
□ 



Following the approach in Part [X] we introduce an e-nct 
Af e in the isometries of V and obtain up to order e 



L(p)=2 E <r(V)|det 

VeJV e 



where we defined 

<r<y)= E l-yix. 

i\...i N 

Consider now the sum 



■ »jv| 2( 5p' 



E = E toi, 

WGAf, ii—i N 



(C2) 



(C3) 



(C4) 



Since it takes a finite value, for any r 3> e there must exist 
an r-region 7£ r C 7V e containing V" such that V V G 7Z r 
the relation II V - V"|| < r holds with 



E ff W 



> 



Vv. 



(C5) 



where and Vjv e are the number of points of the re- 
spective regions in the e-grid. For sufficiently small r we 
consider only the following part of the sum in (C2) 



2 J2 er(V)|det( 

V6K r 

= 2 £ (t(V)| det^" 1 ^)! (1 + 0(0) . (C6) 



If we set Q = V", this is lower bounded by 2(1 
0{r))VnjV Me , i.e., 



L(p)> 2(l + 0(r)) 



(C7) 



which is small but positive for sufficiently small r. □ 



Appendix D: Matrices with LRLE for D = 3 and 

dim(S) = 6 

For D = 3 in the main text it has been shown that the 
case of dim(<S) = 3 corresponds to a trivial extension of 
the matrices for D — 2, Af =2 = (j-jE/*, to D = 3. More- 
over, it has been noted that larger possible dimensions 



of S are only 6 and 8, the latter implying AiA\ oc I. 
For the case of dim(<S) = 6 it has been shown that 
S = span{V", W, W\ V , W, W% where V = \x)(x\ - 
|y)(y|) W = \x)(y \ and analogous definitions with \x') and 
\y'). The reason for the introduction of V and W of the 
same type as V and W is that S is closed under any £j, 
which maps, e.g., V to V% = A\ \ x) (x \ Ai — A\\ y)(y For 
some Si Vi must happen not to be in span{V, VF^, 14^^} 
if dim (5) > 3. We see that in this case A\\x) and A\\y) 
are new orthogonal vectors with the same norm (due 
to tr(Vi) = tr(Wj) = 0). From them we obtain the 
orthonormal vectors \x') = [(x^iA^x)] -1 / 2 A\\x) and 
\y') = [(ylAiAlly^-^Ally). Thus, A t applied from the 
right maps all vectors contained in r = span{|x), \y)} 
to vectors in r' = spanjja;'), \y')}, while preserving 
their relative lengths, or equivalently, the angles between 
them. 

As noted in the main text, the existence of another 
r" ^ t, t' (supporting V" , W" and W"t) is excluded for 
dim(S) = 6. Consequently, any Ai applied from the right 
has to map r to r or r' and r' to r or r' in such a way that 
all angles between vectors lying in one of those subspaces 
are preserved, r and t' ^ r at this point can be arbi- 
trary two-dimensional linear subspaces of C 3 , whereas we 
assume their intersection to be along \y) = \y r ). The jus- 
tification for this choice is that for any orthonormal pair 
la:'), W) of vectors in r' span{V, W, W'^} is the full 
subspace of traceless matrices contained in r' x r'. After 
choosing w.l.o.g. |a;) := |1) and \y) := |2), we can make 
the ansatz 



which results in bi = V i = 1, d. Now, the matrix 



1 
1, =7, ( 1 | I';. 

ai bi c 



(Dl) 



7i, a i} bi, a € C and Ui is a unitary. We see that indeed 
tr(Aty^) = 0, i.e., (1|A,4|1) = (2|A,A+|2) = | 7l | 2 and 
ti(A\WAi) = 0, i.e., (21^411) = arc fulfilled. The 
normal vector of r is \n) = |3), and the one of r' is of 
the form \n') = r|l) + s|3), since r and r' intersect in the 
|2)-axis. Because of ti(A\SAi) = V S £ S, it follows 
that A t A\ e S = span{I, |n)(n'|, |n')(n|}. Therefore, 



A 4 A\ = a i I + i (\3)[r{l\ + 8( 



+ [r*|l) + **|3)](3|), 
(D2) 



Mi 




(D3) 



applied from the right maps r to r and r' to, say, f pre- 
serving the angles between vectors lying in one of them. 
As a result, there are four possible cases for the action of 

Ui 



a) Ui 

b) Ui 

c) Ui 

d) Ui 



(t, f) (r, r), 
(r, f) -+ (r, r'), 
(r, f) (/, r), 
(r, f) (/, r'). 



(D4) 
(D5) 
(D6) 
(D7) 



a) and d) imply f = r and hence a, = e JV and Cj = 0, 
whereas b) and c) mean f ^ t. Therefore, in the latter 
case the intersection of t and f is the |2)-axis. Since it 
has to be mapped by Ui to the intersection of r and r', 
which is likewise the |2)-axis, |2) must be an eigenvector 
of Ui. Then, it follows that Ui has a unitary action in 
the |l)-|3)-plane. Hence, b) corresponds to f = r', and 
therefore dj =0, c, = e 4< ^. Last, in c) double application 

of would map f r ^> r', wherefore Mi has to 
carry out a reflection of r' on the r-plane, i.e., ttj = 
and Cj = — 1. We thus obtain either cij = and = e lv 
or ai = e lv and Cj = V i = 1, d and conclude that 
r'=span{|2),|3)}. 

The requirement that Ai maps r to t or r' and r' to 
t or r' also holds for A 2 , from which we deduce that 
Ui must contain only one non-vanishing entry per row 
and column (with magnitude 1). All results combined 
together yield 



Ai oc [e^\l) + e^\3)]{l\ + e 1 ^ \2)(m\ 



with l,m=l, 2, 3, I ^ m. 



(D8) 
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